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Abstract 

A  parallel  development  of  the  semlclasslcal  and  quantum  statistics  of 
multi -spatlotemporal  mode  direct,  homodyne,  and  heterodyne  detection  using  an 
Ideal  (except  for  Its  sub-unity  quantum  efficiency)  photon  detector  Is  presented. 
Particular  emphasis  Is  placed  on  the  latter  two  coherent  detection  configurations. 
The  primary  Intent  is  to  delineate  the  semlclasslcal  theory's  regime  of  validity 
and  to  show,  within  this  regime  of  validity,  how  the  quantum  theory's  signal 
quantum  noise,  local  oscillator  quantum  noise,  the  quantum  noise  incurred 
because  of  sub-unity  detector  quantum  efficiency,  plus  (for  heterodyning  only) 
image  band  quantum  noise  produce  the  quantitative  equivalent  of  the  semlclasslcal 
theory's  local  oscillator  shot  noise.  The  effects  of  classical  fluctuations 
on  the  local  oscillator,  and  the  recently  suggested  dual-detector  arrangement 
for  suppressing  these  fluctuations,  are  treated.  It  is  shown  that  previous 
studies  of  this  arrangement  have  neglected  a  potentially  significant  noise 
contribution. 


2 


f 


I.  INTRODUCTION 

In  coherent  optical  detection  [1]  •  [3],  the  optical  field  to  be  measured 
Is  combined  on  the  surface  of  a  photodetector  with  the  field  of  a  strong 
local  oscillator  laser  whose  center  frequency  Is  offset  by  an  amount  Av  from 
that  of  the  signal  field.  The  detection  scheme  Is  referred  to  as  optical 
homodynlng  If  dv»0,  and  optical  heterodyning  If  dv  "vIfr>0,  with  vIF  being  the 
Intermediate  frequency  In  the  latter  case.  For  both  schemes,  electrical 
filtering  of  the  photocurrent  Is  used  to  select  the  beat  frequency  components 
In  the  vicinity  of  dv,  yielding  an  output  that  contains  a  frequency  translated 
replica  of  the  signal  field  components  that  were  coherent  In  space  and  time 
with  the  local  oscillator  field.  Heterodyne  detection  Is  now  widely  employed 
In  coherent  CO^  laser  radars  [4],  [S],  and  Is  being  vigorously  researched  for  use 
with  semiconductor  Injection  lasers  In  fiber  optics  [6]  -  [8]  and  space 
communications  [9],  [10].  Performance  analyses  In  these  areas  routinely  employ 
the  semlclasslcal  statistical  modal  for  photodetection,  which  Implies  that  the 
fluctuations  observed  In  coherent  optical  detection  with  signal  and  local 
oscillator  fields  of  perfect  amplitude  and  frequency  stability  comprise  an  additive 
white  Gaussian  noise,  representing  local  oscillator  shot  noise. 

It  has  long  been  known  [11]  that  the  semlclasslcal  statistics  for  photodetec¬ 
tion  are  quantum  mechanically  correct  only  when  the  total  field  illuminating  the 
detector  Is  In  a  Glauber  coherent  state  or  a  classically  random  mixture  of  such 
states.  Inasmuch  as  ordinary  light  sources,  Including  lasers  and  light  emitting 
diodes,  obey  this  classical  state  condition,  there  Is  no  need  to  abandon  the 
semlclasslcal  approach  In  the  vast  majority  of  photodetection  sensitivity 
calculations.  However,  non-classical  light  has  been  generated  via  resonance 
fluorescence,  as  confirmed  by  observations  of  its  photon  anti-bunching  [12]  and 
sub-Polssonian  behavior  [13]  In  direct  detaction.  Moreover,  there  is  great 
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theoretlcal  Interest  In  squeezed  states  (also  called  two-photon  coherent  states) 

[14],  [IS],  which  are  non-classical  states  of  considerable  potential  for  optical 
communications  [16]  -  [19]  and  precision  measurements  [20],-  [23],  For  these 
states,  the  quantum  theory  of  photodetection  Is  essential,  and  coherent  optical 
detection  schemes  are  the  most  Interesting.  / 

In  [18],  Yuen  and  Shapiro  developed  the  quantum  descriptions  of  single- 

/ 

detector  optical  homodyne  and  heterodyne  receivers'.  They  employed  a  quasimono- 

chromatic  approximation,  and  assumed  a  coha rent/state  local  oscillator, 

/ 

corresponding  to  perfect  local  oscillator  amplitude  and  frequency  stability. 

Within  these  limitations,  complete  statistics  for  multi -spatlotemporal  mode  detection 
are  available  from  [18].  More  recently  [2«J],  the  fact  that  photodetectors  respond  ■ 
to  photon  flux  rather  than  power  [24]  hai  been  used  to  relax  somewhat  the 
quaslmonochromatlc  approximation  In  [18].  Thus,  were  high  power,  highly  stable 
local  oscillator  lasers  available  at  all  wavelengths  of  Interest,  the  quantum 
photodetection  theory  of  [18]  would  provide  a  sufficiently  general  foundation  for 
all  optical  homodyne  and  heterodyne  sensitivity  calculations.  Unfortunately, 
such  Is  not  the  case.  1 

□riven  by  heterodyne-detection  problems  arising  from  the  excess  noise  of 

semiconductor  Injection  lasers,  Yuen  anil  Chan  [25]  proposed  a  dual -detector 

\ 

arrangement  for  coherent  optical  detection,  akin  to  the  balanced  mixer  concept  of 
microwave  technology  [26],  [27],  They  gaye  a  direct  quantum  analysis  of  single¬ 
mode  dual-detector  homodynlng,  showing  tha^  local  oscillator  quantum  and  excess 
noises  can  be  balanced  out,  hence  allevlatlr^  Injection  laser  problems  that  would 
have  plagued  a  single-detector  system,  In  subsequent  work  by  Chan  and  his 
collaborators,  the  basic  dual-detector  excess  riolse  cancellation  concept  was 
demonstrated  experimentally  [28],  and  a  variety  \f  non-ideal  device  effects  (quantum 
efficiency  mismatch,  etc.)  were  analyzed  using  sentyclasslcal  multi -temporal  mode 
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techniques  [29].  Also,  Schumaker  [30]  has  shown  that  the  dual-detector  single- 
mode  homodyne  arrangement  is  better  than  single-detector  homodyning  for  making 
non-classical  squeezed  state  observations,  as  a  result  of  Its  ability  to  cancel 
out  local  oscillator  excess  noise. 

Because  the  quantum  treatments  In  [25]  and  [30]  are  confined  to  single-mode 
situations,  and  the  multi-mode  results  In  [28],  [29]  are  In  essence  semlclasslcal , 
there  Is  as  yet  no  fully  quantum  treatment  of  multi -mode  dual -detector  coherent 
optical  reception.  This  paper  will  develop  such  a  model  by  generalizing  the 
results  of  [18].  Simple  explicit  representations  for  all  of  the  relevant  output 
terms  In  coherent  optical  detection  with  a  strong  but  classically  random  local 
oscillator  field  will  be  derived.  It  will  be  seen  that  the  previous  dual¬ 
detector  analyses  [25],  [28]  -  [30]  neglect  excess-noise  modulation  of  the  signal 
and  quantum  noise  terms,  and  the  first  of  these  modulation  effects  may  significantly 
degrade  output  slgnal-to-nolse  ratio  In  some  circumstances.  Moreover,  because 
of  the  calculatlonal  power  afforded  by  [18],  our  rather  general  quantum  results  are 
more  directly  comparable  with  those  of  the  multi-mode  semlclasslcal  theory  than 
are  the  more  limited  results  of  [25],  [30].  Indeed,  that  comparison  Is  the 
primary  purpose  of  this  paper. 

The  paper's  core,  Section  II,  Is  a  parallel  development  of  the  semlclasslcal 
and  quantum  statistics  of  mul tl-spatlotemporal  mode  direct,  homodyne,  and  heterodyne 
detection  using  an  Ideal  (except  for  Its  sub-unity  quantum  efficiency)  photon 
detector.  The  formulation  therein  for  the  coherent  optical  detection  schemes  will 
assume  perfectly  stable  local  oscillators  In  the  semlclasslcal  models,  and  the 
corresponding  coherent  state  local  oscillators  In  the  quantum  models.  We  use 
Section  II  to  delineate  the  semlclasslcal  theory's  regime  of  validity,  and  to  show, 
within  this  regime,  how  the  combination  of  the  quantum  theory's  signal  quantum 
noise,  local  oscillator  quantum  noise,  the  quantum  noise  Incurred  because  of 
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sub-unity  detector  quantum  efficiency,  plus  (for  heterodyning  only)  Image  band 
quantum  noise  produce  the  quantitative  equivalent  of  the  semlclassical 
theory's  local  oscillator  shot  noise.  In  Section  III  we  address  coherent 
optical  detection  with  classically  random  local  oscillators,  The  technique  of 
Iterated  expectation  Is  used  to  readily  obtain  both  semlclassical  and  quantum 
results  for  this  case.  Single-detector  and  dual-detector  systems  are  considered, 
and  our  results  are  compared,  In  the  case  of  dual -detector  quantum  homodynlng  to 
those  of  [25],  [30],  Finally,  In  Section  IV  we  briefly  discuss  the  Implications 
of  our  work  for  squeezed  state  generation  experiments,  which  is  the  application 
that  motivated  our  analysis, 

II.  SEMI CLASS I CAL  VS,  QUANTUM  PHOTOOETECTION 

The  central  element  of  all  the  photodetectlon  configurations  we  will  consider 
Is  shown  In  Fig.  1.  It  Is  a  surface  photoemlttar  with  active  region 

x  •  (x,y)  «  Ad  In  the  z«0  plane,  Illuminated  by  a  quaslmonochronomatlc  (center 
frequency  vQ)  paraxial  scalar  electromagnetic  wave  from  the  half  space  z<0  over 
an  observation  time  Interval  t  c  T.  This  detector  Is  assumed  to  have  a  constant 
quantum  efficiency  n  over  the  frequency  band  containing  the  Illuminating 
field.  The  output  of  the  detector  Is  a  scalar  current  density  J(x,t)  for 
X  C  Ad,  t  c  r.  As  will  be  described  below,  the  field  characterization  we 
must  employ  for  the  Illumination  Is  either  classical  or  quantum  mechanical,  according 
to  whether  semlclassical  or  quantum  photodetectlon  statistics  are  sougnt.  Although 
we  shall  neglect  Internal  time  constant  and  noise  effects,  which  are  present  In 
real  detectors,  our  direct  detection  results  will  be  applicable  to  photomultiplier 
tubes  (for  which  the  current  gain  permits  Internal  noise  to  be  overcome)  at 
post-detection  bandwldths  up  to  the  reciprocal  anode  response  time  of  the  tube, 
Furthermore,  our  results  will  be  applicable  to  coherent  optical  detection  systems 
using  semiconductor  photodiodes  (for  which  the  mixing  gain  overcomes  the  Internal 
noise)  up  to  the  post-detection  bandwidth  of  the  detector.  No  particular  loss  of 


generality  Is  entailed  by  the  use  of  scalar  rather  than  vector  fields,  with 
the  caveat  that  all  the  coherent  optical  detection  work  herein  presumes  that 
the  actual  signal  and  local  oscillator  fields  are  co-polarlzed.  Finally,  by 
appropriate  spatial  Integrations,  we  can  collapse  our  current  density  observation 
to  photocurrent  observations  for  a  single  detector  or  a  multiple-detector 
array. 

A.  Direct  Detection 

In  direct  detection,  the  electromagnetic  field  to  be  measured  comprises  the 
entire  Illumination,  and  the  basic  observation  quantity  Is  the  current  density  J(7,t). 
Semi classical  Model  Let  E^(7,t)  be  the  positive-frequency  complex 

field  (V/m  units)  associated  with  the  classical  scalar  electric  field  Incident 
on  the  detector,  l.e.,  E^(7,t)  Is  the  analytic  signal  of  this  electric  field. 
Because  of  our  quaslmonochromatlc  assumption,,  the  Fourier  transform  1  of  E^ 

E^(7,v)  -  jdt  E(+*(7,t)  (1) 

Is  non-zero  only  for  |v  -  vfl|  <.B,  where  the  bandwidth  B  is  much  less  than  the 
center  frequency  vQ.  Because  of  our  paraxial  assumption,  the  short  time  average 
power  density  falling  on  the  point  7  at  time  t  Is 

1(7, t)  -  (ce0/2)  E(_)(7,t)E(+)(7,t),  (2) 

where  c  1  r.  the  speed  of  light,  eQ  the  permittivity  of  free  space,  and  F-")  s(E^)* 

Is  the  negative- frequency  complex  field,  with  denoting  complex  conjugate. 

The  standard  semlclasslcal  photodetection  model  [31],  In  our  notation, 
presumes  that  0(7, t)  Is  a  conditional  space-time  Poisson  Impulse  train  with 
rata  function  u(7,t)  ■  RI(7,t)/e  where  e  Is  the  electron  charge,  and  P. 

Is  the  detector's  responslvlty  (A/W  units)  at  the  Illumination's  center 
frequency  vQ  .  This  means  that: 
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1)  the  current  density,  which  is  of  the  form 

J(x.t)  *  ^ei(x-xn)6(t'tn) ,  (3) 

n 

has  shot  effect  noise,  l.e.,  it  consists  of  instantaneous 
emissions  of  an  electron  charge  e  at  the  random  space-time 
points  ((7n,tn)  :  ”neAd,  tneT>; 

2)  conditioned  on  knowledge  of  the  rate  function  {u(x,t)  :  xeAd,teT}, 
the  number  of  photoemissions  occurring  within  a  spatial  region 
A'CAd  during  a  time  Interval  7“  C  T  Is  a  Poisson  random  variable 

with  mean  value  dx"  dt  u(7,t)i 
A'  T* 

3)  conditioned  on  knowledge  of  (uCx.t)  :  7tAd  ,  teT>,  the  photoemissions 
occurring  In  disjoint  spaial  regions,  A',  A"CAd  are  statistically 
independent  processes, 

4)  conditioned  on  knowledge  of  Cu(x.t):  7*Ad,tcT),  the  photoemissions 
occurring  In  disjoint  time  Intervals  7'  ,T"CT  are  statistically 
Independent  processes. 

Even  though  the  semlclasslcal  theory  of  photodetection  employs  classical 
fields,  It  Is  customary  to  recognize  In  this  theory  that  light  of  frequency  v 
Is  quantized  into  photons  of  energy  hv  ,  where  h  Is  Planck's  constant.  Thus, 
for  the  quaslmonochromatlc  case  at  hand,  the  responslvlty  is  ordinarily  written 
as  R  ■  en/h\io,  In  terms  of  the  detector's  quantum  efficiency  n  and  the 
photon  energy  at  the  field's  center  frequency,  so  that  u(7,t)  ■  nl(7,t)/hv  . 
In  fact,  because  we  are  concerned  with  detectors  that,  quantum  mechanically, 
respond  to  photon-flux  density  rather  than  power  density  [24],  [22],  it  Is  more 
proper  to  write 

u(7,t)  «  nlph(x,t), 


(4) 
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where  I  ^  Is  the  classical  photon- flux  density 


Iph(x,t)  ■  E#(x\t)E(x,t) 


(5) 


1/2 

obtained  from  the  (photons/s)  /m  units  positive-frequency  complex  field 


E(x,t)  -  dv(c«„/2hv)  1/2  S(+)(x,v)e‘j2,rvt, 


(6) 


For  all  practical  purposes  in  the  semi  classical  theory,  with  quaslmonochromatlc 
light  we  can  use  u(7,t)  ■  n I (x-* t)/hvjQ  and  y(x,t)  ■  nl  h(x,t)  Inter¬ 
changeably.  This  amounts  to  using  v  :  In  the  square-root  term  of  (6), 

an  approximation  whose  validity  Is  guaranteed  by  (1).  In  the  quantum  theory, 
even  with  quaslmonochromatlc  light,  It  Is  critical  to  employ  the  photon-flux 
density  formulation,  see  [22]. 

Quantum  Model  In  the  quantum  photodetectlon  theory,  the  classical  positive- 

frequency  complex  field  E^(x,t)  Is  replaced  by  a  positive-frequency  field 

operator'  E'  ;(x,t),  whose  quantum  state  Is  specified  by  a  density  operator  o. 

The  quaslmonochromatlc  and  paraxial  conditions  of  the  semi  cl  ass  leal  theory 

become  conditions  on  the  density  operator,  namely,  that  the  excited  (non-vacuum 

state)  inodes  of  E^(x,t)  He  at  frequencies  within  B  of  \iQ  and  propagate 

at  small  angles  to  the  z  axis.  As  In  [17],  [18],  [22],  we  shall  regard  the 

current  density  J(7,t)  as  a  classical  quantity,  corresponding  to  the  macroscopic 
2 

output  of  the  quantum  measurement  performed  by  the  detector  of  Fig.  1  on 
the  field  V  ;(x,t).  To  provide  an  explicit  representation  of  this  quantum 
measurement,  we  must  first  develop  the  auantum  effective  photon-flux  density. 

Let  us  convert  ;(x,t)  to  a  photon-units  field  operator  by  defining 
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(Cf.  Eqs.  (1),  (6)  ) 


iw(7.v) 


dt  E^(x\t)  e 


J2irvt 


and 


E(x,t) 


dv  (c*  /2hv)1/2  E^(x,v)e"j2lTvt 


(7) 

(8) 


Equation  (8)  defines  the  same  basic  field  operator  used  In  [16]  -  [18],  [22]. 

Me  can  made  a  modal  expansion 

E(x,t)  -  5nCn(x.-t).  xeAd,  teT  (9) 

of  this  operator,  where  Can>  are  modal  annihilation  operators  satisfying  the 
commutation  rules 

tVSJ  CvO**™- 

and  (cn)  are  a  complete  orthonormal  set  of  classical  functions  over 
xeA^,  tef.  In  Eq.  (10),  the  tan+)  are  the  adjolnts  of  the  { an } ;  they 
are  modal  creation  operators. 

For  a  detector  of  sub-unity  quantum  efficiency  we  must  adjoin  to  (9)  a  fictitious 
field 


Eyj^(x,t)  *  l  C^5p(x,t) ,  XeA^,  tet , 


(ID 


where 


“n> 


n 

are  modal  annihilation  operators  that  commute  with  { an } 
viz. 


and 
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The  fields  E  and  Eyac  are  quantum-mechanlcally  Independent,  with  the  latter 

A  A 

having  all  Its  modes  In  the  vacuum  state.  In  terms  of  E  and  Eyac 
the  effective  photon-flux  density  operator  for  the  detector  Is 

i'h(x,t)  •  E1+(7,t)E'(7.t)  ,  (13) 

with 

E'(x.t)  3  n1/2E(7,t)  +  (l-n)1/2Evac(7,t)  .  (14) 

The  representation  theorem  of  quantum  photodetectlon  [18,  theorem  1]  can 
now  be  stated  (In  our  notation)  as  follows.  The  classical  current  density 
J(7,t)  obtained  from  photoemlssl ve  detection  measures  the  quantum  operator 

3(7, t)  ■  ai'h(7,t)  ,  (15) 

l.e.,  It  Is  proportional  to  the  effective  photon-flux  density.  In  somewhat 
more  detail  this  means  (cf.  the  semlclasslcal  case): 

1)  the  current  density  obeys 

J(7,t)  ■  l  e  6(7  -7„Wt-tJ  ,  (16) 

“  n  n 

n 

so  It  Is  still  a  collection  of  Instantaneous  emissions  of  an  electron 
charge  at  random  space-time  points  C(7  ,tn)}  ; 
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2)  if  F(j(x,t))  is  an  arbitrary  functional  of  the  current  density,  then 
the  classical  average  of  this  random  variable  <F(J(x,t))>  equals 
the  quantum  average  tr(p 'F(J(x,t))),  where  tr  denotes  trace  and 

o'  ■e®Pvac  gives  the  joint  density  operator  for  E  and  Evac  In 
terms  of  the  density  operator  g  for  E  and  the  vacuum-state 

pvac  “  ®  10-01  density  operator  for  Evac 

Note  that  we  cannot  dispense  with  the  vacuum  state  field  Eyac  unless  n«1, 

even  though  Its  average  value  obeys  tr(Pw,„  E„„.Cx.t))  ■  0  regardless  of  the 
value  of  n  •  This  Is  because  the  zero-point  fluctuations  (vacuum-state  quantum 
noise)  In  Evac  can  contribute  to  F(J(x,t)).  Indeed  the  noise  In  J(7,t)  has 
nothing  to  do  with  the  shot  effect  associated  with  the  discreteness  of  the 
electron  charge,  Rather,  it  Is  the  quantum  noise  in  £'  being  observed  through 
measurement  of  the  effective  photon-flux  density. 

Comparison  Let  us  suppose  that  the  density  operator  for  E  is  a  classical  state3, 

l.a. , 


2  * 
da  P(a;a  )  |a»<a| 


(17) 


with 

*  an (18) 

defining  the  multi-mode  Glauber  coherent  states  of  the  field  E  in  terms  of 
the  modal  expansion  (9),  and  P(a;a  )  being  a  classical  probability  density 
function  P(aia  )_>0,  [d^a  P(a;a  )  ■  1,  It  was  shown  in  [18]  that  this  Is  a 

necessary  and  sufficient  condition  for  the  semlclasslcal  statistics  to  be 
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quantitatively  correct,  In  particular,  under  this  condition  the  quantum  model 
predicts  that  J(x,t)  Is  a  conditional  space-time  Poisson  Impulse  train  with 
conditional  rate  function 

u(x.t)  •  o£  (x,t)E(x,t)  (19) 


for 

E(x,t)  a  l  an5n(x,t)  »  <a[E(x,t) |a>  ,  (20) 

n  ■ 

the  average  field  illuminating  the  detector  given  the  state  of  E  Is  the  multi-mode 
coherent  state  |ct  >  , 

To  Illustrate  the  above  behavior,  let  us  examine  the  statistics  of  the 

4 

observed  photon  count 


N 


dt  J(x,t) 
>T 


(21) 


assuming  single-mode  Illumination,  and  n  ■  1.  In  the  semiclasslcal  theory  we 
shall  taka 


E(x.t)  -  Ca/(AdT)1/2]e'j2Tvot.  xeAd,tcT  , 


(22) 


where  a  Is  a  complex-valued  random  variable  with  probability  density  function 
p(a),  Ad  Is  the  area  of  Ad,  and  T  Is  the  duration  of  T.  We  then  obtain 
Handel's  rule  [32]  for  the  probability  distribution  of  N 


d2* 


p(a)  ( 


1 2n 


/n ; )  exp  (- ict  1 2) 


Pr[N  ■  n] 


(23) 
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viz.  N  is  a  conditionally  Poisson  random  variable.  Equation  (23)  gives 
the  mean  and  variance  of  the  observed  photocount  to  be 

;  <N>  ■  d2a  p (a)  |o | 2  •  <|a|2>  (24) 

. 

and 

var(N)  *  <N>  +•  var  (|a|2),  (25) 

respectively,  where  the  first  term  on  the  right  In  (25)  represents  shot 
noise  and  the  second  term  on  the  right  In  (25)  represents  excess  noise. 

In  the  quantum  theory  we  let  $,(7,t)  ■  (AdT)"1/,2e"^2ffvot  be  the  only 
excited  mode  In  (9),  so  that  the  density  operator  for  E  is 

p  ■  p.|  ®  1 0><Q |  (26) 

for  p1  the  density  operator  of  mode  1.  We  then  find  for  the  probability 
distribution  of  N  [11],  [17] 

Pr[N  »  n]  ■  <n|p^  |n>  ,  (27) 

where 

a^+a^ |n>  *  n|n>  (28) 

defines  the  photon  number  states  of  the  first  mode  of  E.  If  is  the 
classical  state. 

*  p 

•  da  p(a)  |a><a|  ,  (29) 

with  p(a)  being  the  probability  density  from  the  semlclasslcal  theory 
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(cf.  Eq.  (17)),  then  (27)  reduces  to  (23)  as  expected.  Thus,  In  this  cr^se  the 
semi  classical  theory  Is  quantitatively  correct  In  Its  prediction  of  the  photon 
counting  probability  distribution.  It  Is  nevertheless  physically  Incorrect,  In 
that  It  ascribes  the  photon  counting  fluctuations  to  shot  noise,  whereas  they  are 
actually  a  manifestation  of  the  Illumination  field's  quantum  noise.  For 
example,  were  p1  *  |k><k|  where  |k>  Is  the  k-photon  number  state  (a 
non-classical  state),  then  we  would  get 

Pr[N«n]  ■  Snk  (30) 

from  (27),  whence 

<N>  «  k  ,  (31) 

and 

var(N)  ■  0  ,  (32) 

for  the  photon  count  mean  and  variance.  Here  the  field  state  Is  an  elgenket  of 
our  observation  operator,  so  there  Is  no  uncertainty  In  the  measurement 
outcome.  This  sub-Pol ssonl an  behavior  cannot  be  obtained  from  the  semi  classical 
theory,  because  for  all  p(a)  the  excess  noise  term  In  (25)  will  be  non-negative, 
forcing  var(N)^  <N>  to  prevail. 

B.  Homodyne  Detection 

The  configuration  we  shall  consider  for  single-detector  multi -spatlotemporal 
mode  homodyne  detection  Is  shown  In  Fig.  2.  The  signal  field  to  be  detected  Is 
combined,  through  a  lossless  beam  splitter  of  Intensity  transmission  e,  with  a 
perfectly  stable  local  oscl 11 ato»*  field  on  the  surface  of  the  Fig.  1  photodetector. 
The  resulting  currant  density,  Jhom(x,t),  Is  our  homodyne  detection  output, 
whose  statistics  we  shall  characterize  below.  By  spatial  Integration  of  our  results 
over  the  detector's  active  region  A^,  we  can  use  our  model  to  describe  single- 
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detector  homodynlng;  the  extension  to  dual-detector  homodynlng  will  be  made  In 
Section  III. 

Semi  classical  Model  The  total  classical  photon-units  complex  field  Incident  on  the 
photodetector  Is  given  by5 

E(x,t)  «  c1/2Es(x,t)  +  (l-e)1/2EL0(x,t)  ,  (33; 

for  7sAd  ,  tef,  In  terms  of  a  (potentially  random)  weak  signal  field 

Es(7,t),  and  a  deterministic  strong  local  oscillator  field  ELg(x,t).  The 
latter  has  a  classical  photon-flux  density 

IphLO^)  “  £L0*(x,t)EL0(x,t)  ,  (34) 

that  greatly  exceeds  that  of  the  former 

IphS(x,t)  “  •  (35) 

for  xcAd,tcT  .  Thus,  the  rate  function  driving  the  photodetector  Is,  from 
(4),  (5),  (33)-(35) ,  approximately 

u(x,t)  -  n[(l-c)IphL0(x,t)+2Cc(l-c)]1/2Re(Es(x,t)EL0*(x,t))]  .  (36) 

It  then  follows,  from  the  Central  Limit  Theorem  for  high  density  shot  noise  [33], 
[34],  that  at  very  large  values  of  the  local  oscillator  classical  photon  number 


NphlO 


rdt 


(37) 


the  homodyne  detection  current  density  J(10m(x,t)  Is  a  conditional  Sausslan 
process.  Specifically,  conditioned  on  knowledge  of  the  signal  field 
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vEj(x,t)  :  xcAd,  tcT)  ,  t^ie  sum  t^iree  current  densities: 

1)  a  homodyne-mlxlng  currant  density  signal  term  2en[e(l-c)]1/,2Re(Es(x,t)ELQ*(x,t)) 

2)  a  dlmct-detectlon  local  oscillator  bias  current  density  en(l-c)IphLg(7,t)  ; 
and 

3)  a  local  oscillator  shot  noise  current  density,  which  Is  a  zero-mean 

spatlotemporal  non-statlonary  white  Gaussian  noise  process  with 

covariance  function 

<Jihot^l,tl^Jshot^2,t2^>  " 

eZn(l-e)lphL0(x‘1  ,t1)i(x1-x2)«(t1-t2)  .  (38) 

In  order  to  connect  the  preceding  multl-spatlotemporal  mode  formulation  with 
more  familiar  single-detector  multi-temporal  mode  results,  let  us  consider  the 
statistics  of  the  single-detector  homodyne  photocurrent 

W‘>  •  j.  d*  W*'4'  •  <») 

*d 

assuming  that 

EL0(x,t)  •  (PL0/^oAd)1/2  e^V  ,  (40) 

corresponding  to  a  normally-lncldent  plana  wave  local  oscillator  of  power  P^g. 

Hera  we  find  that,  conditioned  on  knowledge  of  the  signal  field,  1^  (t) 
comprises  a  signal  current 

1s1g(t)  ■  2en[PLOc(l-e)/hv0Ad]1/2Re(j  dx  E$(7,t)eJ  2’V)  ,  (41) 

Ad 


plus  a  zero-frequency  bias  current 


W  ■  •'><’-« >PL0/hvt  .  («) 

plus  a  zero-mean  stationary  white  Gaussian  process  shot-noise  current 
1Sh0t(t)  with  spectral  density  e1b1a$  (A2/Hz)  , 

The  signal  current  Is  a  frequency-translated  (to  baseband)  replica  of  the 
normally-lncldent  plane  wave  component  of  Eg(7,t)  that  Is  In  phase  with  the 
local  oscillator  field.  The  bias  current  Is  the  zero-frequency  photocurrent 
produced  by  the  local  oscillator  field.  The  noise  current  Is  the  local -oscillator 
shot  noise,  whose  spectrum  follows  the  well  known  Schottky  formula  [35]. 

Quantum  Model  In  the  quantum  model,  £q.  (33)  becomes  an  operator-valued  expression 

E(7,t)  -  t1/2Es(x,t)  t  (l-e)1/2EL0(x,t),  (43) 

giving  the  field  operator  E  that  drives  the  detector  In  terms  of  the  signal 
field  operator  E$  and  the  local  oscillator  field  operator  ELQ  .  The  density 

A 

operator  g  for  E  Is  assumed  to  be 

g“  Dj  ®  o^q  (44) 

where  oc  Is  an  arbitrary  signal  field  density  operator  and  o,n  *  |a  ><a 

-L0  -L0 

Is  a  multi-mode  coherent  state  local  oscillator  density  operator.  The  latter 
corresponds  to  a  mean  local  oscillator  field 

EL0(7,t)  i  <aLO|ELO(x,t)|aLO>  -  Jj  €n(x,t)  ,  (45) 

when  ^LO  fs  expanded  using  the  mode  set  Un>  as  was  done  for  E  In  Eq.  (9). 


r 
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The  strong  local  oscillator  condition  of  the  quantum  theory, 


tr<psEs+(x\t)Es(7,t))  «  |EL0(7,t) 


Is  assumed  to  prevail  (cf,  Eqs,  (34),  (35)),  with  a  very  large  average  local 
oscillator  photon  number 

Nlq  i  dx  dt|EL0(x,t)'|2»  1,  (4?) 

A  T 

(cf.  Eq.  (37)). 

To  obtain  the  effective  photon-flux  density  operator  measured  by  the  detector 
we  adjoin  to  E  from  (43)  a  quantum-mechanical ly  Independent  vacuum-state 
field  operator  Evac  ,see  Eqs.  (11)  -  (14).  We  can  now  give  a  fully  quantum 
characterization  of  the  classical  homodyne  current  density  JhQm(“,t)  ,by 
translating  the  results  of  [18,  theorem  2]  Into  our  notation.  The  strong  local 
oscillator  condition  Implies  that  this  classical  current  density  measures  the 
quantum  operator 

Jhom(7,t)  -  en(l-c)Ej0(x,t)EL0(7,t)  +  2e[n(l-c )]1/2Re  ([(nc)1/2Es(x,t) 

+  n-n)1/2*Evac(7,t)]  •  6L0*(x,t)  >  .  (48) 

Moreover,  because  nlo»i,  the  local  oscillator  direct  detection  term  In  (48), 

an(l-c)ELg(7,t)ELg(7,t)  ,  yields  classical  observation  values  comprising  a  bias 
current  density  en(1-c ) |EL0(x,t) |  plus  a  local  oscillator  quantum  noise  current 
density,  which  Is  a  zero-mean  spatlotemporal  non-statlonary  white  Gaussian  noise 
process  J|^q(x,t)  with  covariance  function 

<JL0q^xl  ,1:1  *t2>>  " 

[en(l-c)  ]2  [  E'lq (T'i  ,t1  )(2'5(71-x2).i(t1-t2) 


(49) 
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Furthermore,  under  this  same  condition6,  the  second  term  on  the  right  In  (48) 
simplifies  to  a  homodyne-mlxlng  signal  operator  2en[e(l-c)]  '  Re(Ej(x,t)ELg(x,t)) 

plus  a  sub-unity  quantum  efficiency  (n<l)  quantum-noise  current  density.  The  latter 
current  density  Is  a  zero-mean  spatlotemporal  non-statlonary  white  Gaussian 
noise  process  Jvac(><it)  with  covariance  function 

<(W*rV'Wx2,t2)>  “ 

e2n(1-n)(l-e)|6LQ(x1,t1)|2i(x1-x2)«(t1-t2)  ;  (50) 

Jvac  s  Independent  of  JLQq.  Thus,  the  classical  homodyne  current 

density  Jhcm(x,t)  measures  the  operator 

3hom(x,t)  -  en(1-c)|EL0(x,t)|2  +  JLQq(x,t) 

+  2enCc(l-c)]1/Z  Re(Es(x,t)ELJ(x,t))  +  Ovtc(x,t).  (51) 

The  first  term  on  the  right  In  (ST)  Is  the  local  oscillator  bias  current  density, 
the  second  term  Is  the  classical  representation  of  the  local  oscillator's  quantum 
noise  contributed  by  the  en< 1 -c E^q  measurement,  and  the  last  term  Is  the 
classical  representation  of  the  n<l  quantum  noise  contributed  by  the 
2eCn(l-n)(l-c)]1'/2Re(Evac(lT,t)el_Q(x,t))  measurement.  The  signal  field 
contribution  to  cannot  be  simplified  further  without  knowledge  of 

the  density  operator  o$  .  In  general,  this  term  will  contribute  signal  field 
quantum  noise  to  the  homodyne  observation,  as  will  be  seen  below. 

Comparison  To  facilitate  comparison  of  the  semi  classical  and  the  quantum  theories 
of  homodynlng,  we  shall  restrict  our  consideration  to  the  single  detector  case. 

First,  we  need  the  quantum  characterization  of  the  homodyne  photocurrent  (39), 
which  can  be  obtained  by  spatial  integration  of  the  results  just  presented.  We 
assume  a  normally-lncl dent  plane  wave  mean  local  oscillator  field 
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(52) 


EL0(x.t)  -  (PL0/hvoAd)1/2e-j2"V  , 
and  we  find  the  1^  (t)  measures  the  operator 

^hom(t)  "  ^ bias  +  ^Oq^ 
f  2*nCPL0cO-c)/hvoAd]1/2RB([  dx!s(x,t)  « J2™ot)  +  1yac(t)  _  (S3) 

i 

Here,  Is  given  by  (42),  and  1^  and  1yac  ar«  statistically 

Indtpendant  zero-mean  stationary  white  Gaussian  noise  processes  with  spectral 
densities  «n(l-e)1b1as  and  «(1*n)1bia|  .respectively. 

Physically,  1b1as  Is  the  local  oscillator  bias  current,  1^  is  the  local 
oscillator  quantum-noise  current,  and  1yac  Is  the  n<l  quantum-noise  current. 
Equation  (53)  differs  from  the  semlclasslcal  description  In  two  respects:  the 
homodyne-mixing  signal  tertn  Involves  the  quantum  field  operator  Es  rather 
than  the  classical  field  Eg  ;  the  noise  In  the  homodyne  observation  Is  a 
combination. of  local  oscillator  quantum  noise,  n<l  quantum  noise,  and  signal 
quantum  noise,  rather  than  simply  being  local  oscillator  shot  noise.  We  know, 
from  the dl rect  detection  discussion,  that  the  semlclasslcal  photodetection  model 
Is  quantitatively  correct  If  the  density  operator  p  for  the  field  E  Illuminating 
the  detector  represents  a  classical  state.  This  situation  occurs,  under  (44), 

If  and  only  If  ps  ,the  signal  field  density  operator,  Is  a  classical  state 


for  |as>  the  multi-mode  signal  field  coherent  state  In  modal  expansion  of 
Eg  similar  to  Eq.  (9),  with  Pg  being  a  classical  probability  density.  Whan  (54) 
applies,  the  homodyne-mlxlng  signal  term  In  (53)  can  be  given  a  classical 
representation  akin  to  that  employed  for  the  Eyac  mixing  term  In  going  from 


(48)  to  (51).  In  particular,  for  a  classical  signal  field  state,  the  quantum 
theory  of  homodynlng  predicts  that 

•  ’bias  +  W‘>  + 

SsnCP^eO-O/hv^f^RtC  dx£j(x,t}  t  ’2w;ot)  + 

V*>  *1„e(t)  .  (55) 

where 

Es(x,t)  i  <aj|Ej(>r,t) |«s>  (56) 

Is  the  classical  mean  signal  field  when  the  state  of  Is  |«s>  .  and 
lSq(t)  Is  a  zero-mean  stationary  white  Gaussian  noise  current  of  spectral 
density  ene1b1as  that  Is  statistically  Independent  of  1^  and  1vac 

The  classical  field  Es  Is,  In  general,  a  random  process  with  probability 
density  Ps(as;as*)  In  modal  expansion  form.  The  current  1Sq(t)  Is  the 
classical  representation  of  the  coherent  state  signal  field  quantum  noise  as 
observed  through  the  measurement  operator  (S3).  Note  that 

<'(«)  <  \oqtt)  +  1Sq(t)  *  1v,c<t)  (57) 

Is  a  zero-mean  stationary  white  Gaussian  noise  process  of  spectral  density 
e1b1aj  ,1 n  quantitative  agreement  with  the  semlclasslcal  1jhot(t)  result. 

Of  course,  the  Interpretation  of  the  origin  of  the  noise  In  homodynlng  Is 
different  In  these  two  theories.  Local  oscillator  shot  noise  Is  a  semlclasslcal 
fiction;  the  noise  seen  In  homodyne  detection  (with  an  Ideal  local  oscillator) 
Is  local  oscillator  quantum  noise,  plus  n<1  quantum  noise,  plus  signal  quantum 


noise.  Moreover,  In  the  limit  ,ie  -*■  1  with  n ( 1  -c ) N^q  >>  1  .the  two  former 
contributions  disappear,  and  homodynlng  gives  a  direct  quantum  measurement  of 
the  signal  field  component  that  Is  coherent  In  space  and  In  phase  with  the 
local  oscillator  [18].  It  Is  this  characteristic  that  makes  homodynlng 
attractive  for  squeezed  state  applications  [16]  -  [19]. 

C.  Heterodyne  Detection 

The  configuration  for  single-detector  multl-spatiotemporal  mode  heterodyne 
detection,  shown  In  Fig,  3,  mimics  that  employed  for  homodyne  detection.  The 
only  differences  are  that  the  signal  field  Is  centered  at  frequency  vQ  +  vjp  , 
the  local  oscillator  Is  centered  at  frequency  vQ  ,and  passband  filtering  of 
the  current  density  1$  used  to  select  beat  frequency  components  In  the  vicinity 
of  the  IF  frequency  vIp  (viPT»1  will  be  assumed).  The  bandwidth  B  of 
the  signal  field  will  be  taken  to  be  much  less  than  vjp  ,and  we  shall  concern 
ourselves  with  characterizing  the  statistics  of  the  current  density 
Jhtt(x,t).  The  results  we  need  are  easily  developed  by  Injecting  the  frequency 
offset  vjp  Into  the  preceding  homodyne  work, 

Semlclasslcal  Model  In  Eq.  (33)  let  us  make  the  frequency  offset  of  the  signal 
field  explicit  by  writing 

Es(x,t)  ■  Fs(x,t)e‘J&ff(VvIF)t  ,  (58) 

where  F^  Is  a  baseband  complex  signal  field  of  bandwidth  3  ,  The  results 

following  (37)  now  provide  the  semlclasslcal  statistics  for  heterodyning,  namely, 
conditioned  on  knowledge  of  the  baseband  signal  field  {Fs(7,t)  : 

7cAd,  tcT}  ,  J^#t(7,t)  Is  the  sum  of  three  current  densities: 

1)  a  heterodyne-mixing  current  density  signal  tsrm 

2en[c(l-c)]V2Re(Fs(7,t)e‘J2,T(VJIF)tEL0*(x,t).)  •, 
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2)  a  direct-detection  local  oscillator  bias  current  density 
en(1-c)IphL0(x..t)  }  and 

3)  a  local  oscillator  shot  noise  current  density  charac-t8r-,zetl 

by  (38). 

The  single-detector  heterodyne'  photocurrent 

1het(t)  "  L  d*  Jhet(7<t)  *  (59) 

*d 

assuming  EL0  1s  given  by  (40),  then  comprises  a  signal  current 

1s1g(t)«2en[PL0c(l-.)/  hv0Ad]1/2Re(f  d5frs(x,t)e‘j2,,yIFt)  ,  (60) 

kd 

plus  a  bias  current  1b1as  from  (42),  plus  a  zero-mean  stationary  white  Sausslan 
process  shot-noise  current  1$hot(t)  with  spectral  density  e1b1aj, 

The  heterodyne  current  (59)  'Is  thus  a  frequency  translated  (from 
vo  +V1F  t0  yiF^  vir8,|on  of  honnally  Incident  plane  wave  component  of 
Es  plus  the  usual  bias  and  shot  noise  terms.  Because  of  the  frequency  offset 
VjP  between  the  signal  and  the  local  oscillator  fields,  both  the  in-phase  and 
quadrature  (relative  to  the  local  oscillator)  components  of  the  signal  field 
contribute  to  the  output  observations, 

Quantum  Model  Hare  we  suppose  that  the  only  non-vacuum  state  modes  of  the 
field  operator  Es  lie  within  a  bandwidth  B  of  the  frequency  vQ+vjp. 

However,  because  of  zero-point  fluctuations,  the  quantum  version  of  (58)  Is 

Es(x,t)  -  Fs(x,t)e'jaT{VwIF)t  +  FI(x,t)e‘jar(V''lF)t  ,  (61) 

where  Fs  and  Fj  are  baseband  complex  signal  and  Image  field  operators, 
Physically,  the  Image  band,  being  vIF  Hz  below  the  local  oscillator's  frequency, 
contributes  quantum  noise  to  Jhet  even  when  It  Is  unexcited  [18],  [22], 

We  shall  assume  that  Fs  and  F,  are  quantum-mechanlcally  Independent,  with 
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the  latter  having  all  its  modes  In  the  vacuum  stated  We  now  find,  from 
the  quantum  homodyning  work,  that  measures  the  operator 

3hct(x.t)  "  «nO-c)|ELQ(x,t)  |2  +  0LOc,{x,t) 

+  2enCc(l-e)]1/2Re(Fs(x,t)e'J2Tr(vo+vIF}tHL0*(x,t)) 

+2enCc(l-e)]1/2Re(FI(x,t)e’J2lr(vo‘vIF)tEL0*(7.t)J  +  Jvac(x,t)  ,  (62) 


where  E^q,  JLOq  ,  and  Jyac  are  as  given  in  (SI).  We  can  use  the  vacuum-state 
nature  of  Fj  to  obtain  the  classical  representation 

2en[*(l-«)]1/2Re(FI(x,t)e'J2,r(VwIF)tEL0#(x.t)) 

-  JIq(x,t)/21/2  ,  (63) 

where  Jj  is  a  zero-mean  non-statlonary  white  Gaussian  classical  process 
corresponding  to  the  image-band  quantum  noise,  with  covariance  function 


<Jiq(xi,t1)JIq(x2,t2)>  - 

(en)2e(l-t )  |E|_q(xi  ,t^ )  |  i5(x-|-x2)i(t^-t2) 


(64) 


Thus,  the  quantum  description  of  the  single-detector  heterodyne  photocurrent 
1hat(t)  from  (5$)  i s  that  it  measures  the  operator 


W‘>  •  'ki..  "WO  * 


bias 


‘LOq 


Iq' 


vac' 


+  2enCPLQ£(l-c)/HvqAd]1/2Re(|  dxF$(x  ,t)e“J2wvIFt )  , 

Ad 


(65) 
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where  ^blas  '  \oq  '  and  \ac^  are  aS  1n  and  ^q^  1s  a  zel"° 
mean  white  Gaussian  noise  process  (the  classical  representation  of  Image-band 

quantum  noise)  of  spectral  height  er^c  1  bl as ’ 

* 

Comparison  Suppose  the  density  operator  for  Fs  Is  a  classical  state,  l.e., 

its  density  operator  pP  obeys 

r3 


where  Pc  Is  a  classical  probability  density,  and  |a_>  is  the  multi-mode 
G’auber  coherent  state  for  the  modal  expansion 


P-(x,t)  -  l'i  5n(7,t)eJ2,r(Vv!F)t  ,  (67) 

5  n  Sn  n 

with  l  denoting  summation  over  modes  $n  Uelng  within  bandwidth  B  of 
n 

frequency  Here  we  can  ol5ta1n  a  classical  representation  of  the 

p  term  In  (66)  which  reduces  the  quantum  description  of  the  heterodyne 
¥ 

photocurrent  to 

Wt}  "  1  bl  as  +  1L0q(t)  +  f 

+2enLPLQc(l-c)/hvQAdJV2Re(  ]dx  Fs(x,t)e  •j2lTVIFt) 


+  1Sq(t)/21/2  , 


where  the  total  noise  current,  1LQq  +  Iq 


+  1,„/21/2  «■  S„/2V2  +  1 


vac  ’ 


is  a  zero-mean  white  Gaussian  process  with  spectral  density  el,  .  , 
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In  quantitative,  agreement  with  the  semlclasslcal  theory,  and 

Fs(x,t)  s  <as|Fs(7,t)  |as>  (69) 

A 

Is  the  classical  baseband  signal  field  envelope  associates  with  the 
coherent  state  |as>  .  Note  that  half  of  the  Ej  quantum  noise  entering 

1|iet  comes  throu9h  the  s-|9nal  field  operator  F^  and  the  other  half 

A  0 

comes  through  the  Image  field  operator  Fj. 

III.  EXCESS  NOISE  EFFECTS  ANO  OUAL-OETECTOR  OPERATION 

In  this  section  we  shall  extend  the  results  of  Section  II  for  coherent  optical 
reception  to  include  classical  excess  noise  on  the  local  oscillator  field  and 
dual -detector  operation.  It  Is  convenient  to  begin  with  a  presentation  of 
dual -detector  results  In  the  absence  of  excess  noise. 

A.  Dual -Detector  Coherent  Optical  Reception 

Suppose  the  homodyne/heterodyne  configurations  of  Figs.  2  and  3  are  augmented 
by  the  use  of  another  quantum  efficiency  n  detector  on  the  previously  unused 
output  port  of  their  beam  splitters,  see  Fig.  4.  We  take  the  classical  output 
field  for  this  port  to  be 

E(x,t)  -  -  (l-e)1/2Es(7.t)  +  e1/2EL0(x,t)  (TO) 

In  the  semlclasslcal  model,  and  use  the  corresponding  ooerator-valued 
expression  in  the  quantum  model.  Rather  than  treat  the  full  multi -spatlotemporal 
mode  situation,  we  shall  restrict  our  attention  to  the  photocurrents 
1^(t)  and  1 g ( t)  obtained  by  spatial  Integration  of  the  current  densities 
J.j(x,t)  and  J,(7,t)  produced  by  detectors  1  and  2.  We  shall  assume  a 
perfectly  stable  (l.e,  deterministic)  classical  local  oscillator  field 
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ELQ(x,t)  -  FL0(x,t)e'j2,rvot  (71) 

with  baseband  complex  envelope  F^q  In  the  semi  cl  ass  leal  model,  and  a 
Slauber  coherent  state  quantum  local  oscillator  with  mean  field 

tr(pj.gELQ(x,t))  •  F^Cx.tJe'^V  (72) 

with  baseband  complex  envelope  F10  In  the  quantum  model.  Under  these 
conditions  the  results  of  Section  II  can  be  used  to  show  that  the  following 
statistics  apply. 

Homodyne  Detection  In  homodynlng,  the  signal  field  Is  centered  on  uQ,  so, 
because  of  (71), (72),  It  Is  convenient  to  Introduce  baseband  signal  complex 
envelopes  via 

Es(x,t)  -  Fs(7,t)e"j2*vot  (73) 

and 

EsCx,t)  ■  Fs(7,t)e"^2,rV  ,  (74) 

for  the  semlclasslcal  and  quantum  cases,  respectively.  Now  we  have,  saml- 
classlcally,  that 


1 1 (t)  *  en(l-e) 
*  2an[e(l-t)] 


d7|FLQ(x,t)|2]1/2nshot1(t), 

Ad 


+  e[n(l-e ) 


(75) 
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and 

12(t)  *  enc  dx|FLg(x,t) |2 

Ad  # 

-2enCcO-c)]1/?R®(  *  dx  Fs(x,t)FL0  (x,t)) 

;*d 

+e[nc  dx|FL0(x,t)|2]1/2nsh<Jt2(t)  .  (76) 

Ad 

for  the  homodyne  photocurrents ,  where  nahot1^^  and  n$hot2^  are  statistically 
Independent  Identically  distributed  zero-mean  stationary  white  Gaussian  noise 
processes  of  unity  spectral  density.  Equations  (75)  and  (76)  have  the  usual 
bias  plus  mixing  signal  plus  local -oscillator  shot  noise  Interpretation.  Note 
that  the  beam  splitter  phase  shift  between  the  output  ports  forces  the  mixing 
signals  to  be  180°  out  of  phase.  Also,  the  Independence  of  the  local- 
osclllator  shot  noises  follows  because  they  are  generated  from. deterministic 
Illumination  of  two  different  detectors. 

For  the  quantum  case,  we  have  that  1^(t)  and  12(t)  measure  the  operators 

1^(t)  *  en(l-c)  dx|FLQ(x,t)|2 

Ad  # 

+  2en[c(l-c)]1/2Re(j  dxFs(x,t)FLQ  (x.t)) 

Ad 

+  en(l-e)[j  d7|FL0(x,t)|2]1/2  nLQq(t) 

Ad 

+  e[n(l-n)(l-s)j  dx|FLQ(x.t) i2]1/2nvacl (t)  .  (77) 

Ad 

1 2( t)  *  enc  dx|FLg(x,t) |2 
Ad 

-2en[c(l-e)]1/ZRe(  dx  Fs(x,t)FLQ*(x,t)  ) 

■A  , 


and 
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+ene[  dx|FLQ(x,t) |2]1/2nL(^(t) 

Ad 

+eCn(1-n)if  dx|FL0(x,t) !2]1/2nvac2(  t) .  (78) 

Ad 

where  nLQq(t),  nvacl(t),  nyac2(t)  ar«  statistically  independent  identically 
distributed  zero-mean  stationary  unity- spectrum  Gaussian  processes.  The  familiar 
bias  plus  quantum  mixing  signal  plus  local -oscillator  quantum  noise  plus 
n<l  quantum  noise  interpretation  applies  to  (77),  (78).  As  In  the  semiclasslcal 
model,  the  mixing  term  appears  180°  out  of  phase  In  the  two  photocurrents.  No 
such  phase  shift  appears  on  the  hLfl  term,  as  this  noise  arises  out  of  the 

A 

direct  detection  of  E^q .  Indeed,  except  for  scale  factors,  the  local-oscillator 
quantum  noise  contributions  to  1,(t)  and  12(t)  are  completely  correlated.  The 
n<1  quantum  noises  are,  on  the  other  hand,  statistically  Independent  because  they 
arise  from  different  detectors.  Finally,  when  the  signal  field  is  in  a  classical 
state  these  quantum  results  can  be  shown  to  be  in  quantitative  agreement  with 
the  foregoing  semiclasslcal  formulas. 

Heterodyne  Detection  For  heterodyning  we  use  (58),  rather  than  (73),  to  Introduce 
a  baseband  signal  complex  envelope  for  the  semiclasslcal  analysis.  We  then  find 
that 

i^t)  -  en(l-c)  f  dx|FLQ(7,t)|2 
Ad 

+2en[e(l-«)]1/2Ra(j  dxFs(x  ,t)FL0*Cx  ,t  )a"J2wvIFt:) 

Ad 

+e[n(l-e) j  d7|FLQ(7,t) |231/2n$hotl(t)  , 


(79) 
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12(t)  -  ant  d'x|FL0(x,t)|2 
Ad 

-2an[c(1-c)],/2Ra(j  dxFs(x,t)FLQ*(x,t)e“J2'rvIFt) 

Ad 

+e[ne  dx|FL0(x,t)|2]1/2nshot2(t)l  (80) 

Ad 

with  Intarpret&tlons  as  glvan  following  (75),  (76),  In  the  quantum  case  wa 
usa  (61)  Instaad  of  (74)  and  obtain  tha  measurement  operators 

1-1  (t)  ■  an(1-«) J  dx|FLQ(x,t)  |2 
Md  7 

+2en[t(l-c)]1/2Re(j  dxFs(x,t)F(J[(x,t)e‘J2,TVIFt) 

Ad 

+an(1-c)  dx|FL0(x.t)|2]1/2  nLQq(t) 

d 

+«Cn(l-n)(l-t)j  dx|FL0(x,t)|2]1/2nvacl(t) 

Ad 

+enCc(1-t)2*1  dx|FLQ(x,t)|2]1/2nlq(t)  ,  (81) 

;Ad 

for  detector  1  and 

12(t)  ■  ant  dx|FLQ(x,t) |2 
Ad 

-2anCc(1-t)]I/2Re(j^  dxFjtf.OF^OT.tJa’^IF*) 

+ene[  dx|FL0(x,t) |2]1/2nLQq(t) 

Ad 

+*Cn(T-n)e[  dx| FLQ(x,t) |2]1/2nvac2(t) 

Ad 

-en[e (1-c )2"M  dx|FL0(7,t)|2]1/2nIq(t)  , 

JAd 


(92) 
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for  detector  2.  In  Eqs  (81),  (82)  the  Interpretations  and  comments  following 
(77),  (78)  are  applicable.  The  noise  n^(t),  which  represents  Image-band 
quantum  noise,  Is  another  zero-mean  stationary  unity-spectrum  white  Gaussian 
process.  It  Is  statistically  Independent  of  nLOq(t)  end  nyacj(t)  for 
J  •  1,2,  and  appears  with  a  sign  reversal  In  1^  and  1g  oecause  It  arises 
from  the  mixing  term  Involving  Fj. 

B.  Local -Oscillator  Excess  Noise 

The  extension  of  the  results  of  Section  IIIA  to  Incorporate  classical 
excess  noise  on  the  local  oscillator  Is  extraordinarily  simple,  because  of  the 
form  the  preceding  results  have  been  cast  In.  Specifically,  for  the  semlclassleal 
theory  we  need  only  make  the  baseband  local  oscillator  complex  envelope  flo 
In  (71)  a  comp lex- valued  random  process  with  known  statistics.  Then  the  homodyne 
and  heterodyne  results  of  the  semlclassleal  theory,  npmely  Eqs.  (75),  (76) 
and  Eqs.  (79),  (80),  respectively,  become  conditional  statistics  assuming  FL0 
Is  known.9  Unconditional  statistics  follow,  via  Iterated  expectation  [38], 
from  averaging  over  the  local  oscillator  fluctuations,  as  will  be  Illustrated 
below.  In  a  similar  manner,  classical  local -oscillator  excess  noise  can  be 
Injected  Into  the  quantum  modal  by  making  o^g  a  classical-state  density  operator 
for  which  F^g,  the  average  baseband  local-oscillator  complex  envelope  given  the 
local  oscillator  Is  known  to  be  In  the  multi-mode  coherent  state  |aLg>,  Is.  a 
complex- valued  classical  random  process.  The  quantum  homodyne  and  heterodyne 
results,  Eqs.  (77),  (78),  and  (81),  (82),  respectively,  are  now  conditional 
characterizations  given  F^.  Unconditional  statistics  are  again  obtained  by 
averaging  over  the  local  oscillator  fluctuations.''0 

To  Illustrate  our  excess  noise  results,  and  compare  them  with  relevant  prior 
work  [25],  [28]- [ 30] ,  we  shall  consider  a  single  spatial  mode/mu  1 tl -temporal 
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mode  local  oscillator,  for  which  In  the  semlclasslcal  theory  and  F^g  In 
the  quantum  theory  are  both  of  the  form  (P|_Q(t)/hv0Ad)^2  exp( t ) )  ,  where 
PL0(t)  and  ^g(t)  are  classical  random  power  and  phase  fluctuations.  For 
convenience,  we  shall  assume  that  these  fluctuations  are  the  polar  decomposition 
of  a  stationary  complex-Gauss Ian  random  process.  We  shall  .also  assume  that  the 
signal  field,  In  both  the  semlclasslcal  and  quantum  pictures,  1$  statistically 
Independent  of  the  local  oscillator.  Finally,  we  shall  limit  our  consideration 
to  the  differenced  output  currants  1^(t)  -  12(t). 

Homodyne  Detection  Under  the  preceding  conditions  we  have  the  semlclasslcal 
result 

l^t)  -  12(t)  -  en(l-2c)PLQ(t)/hv0 

+denCc(l-c)PL0(t)/hvoAd]1/2Re(f  dxFs(x,t)e%0(t> ) 

Ad 

^•Cn(l-c)PL0(t)/h,o]1/2nlhotl(t) 

*«CncPLg(t)/hv>0]1/2nJhQt.2(t)  ,  (83) 

and  the  quantum  result 

^(t)  -  12(t)  -  en(l-2e)PL0(t)/hvo 

+4er,re  ( 1  -c )  PL0  ( t )  /hvQ  Ad] 1  /2«e  ( f  dxFs  (7,  t )  aJ*  LO ' ( t : > ) 

Ad 

+en(l-2c)(PLO(t)/hv0)1/2nLOq(t) 

4’e[n(1-n)(l-e)PLg(t)/hv0]1/'2nvacl(t) 

-eCn(l-n)cPL0(t)/hvo]1/2  nyac2(t)  . 


(34) 
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In  both  (83)  and  (84).  the  first  term  on  the  right  equals  a  mean  bias  current  ‘i 

en(1-2e)<Py}(t)>/hv0  plus  a  local  oscillator  power-fluctuation  excess 
noise  en(l-2c)(PLQ(t)-<PL0(t)>)/hvo.  80th  of  these  are  exactly  nulled  when 

1 

the  beam  splitter  Is  50/50,  l.e.,  when  c«  1/2.  The  second  term  on  the  right  In 
(83)  and  (84)  Is  the  homodyne-mlxlng  signal  current*  local  oscillator  randomness 
both  amplitude  and  phase  modulates  this  term.  The  remaining  terms  In  the 
semlclasslcal  result  (83)  are  the  shot  noises,  now  modulated  by  local  oscillator 
power  fluctuations.  The  remaining  terms  In  the  quantum  result  are  the  local 
oscillator  quantum  noise  and  the  n<1  quantum  noises*  these  too  are  modulated 
by  the  local  oscillator  power  fluctuations.  Note  that  when  c  *  1/2  the  local 
oscillator  quantum  noise  contribution  vanishes. 

Let  us  further  specialize  the  quantum  results  by  supposing  that  the  only 

*  *1/2 
excited  mode  of  Fj  Is  the  monochromatic  plane-wave  pulse  ( A^T)  '  for 

xcAd  ,  tcT  ,  and  that  as  Is  the  annihilation  operator  for  this  mode.  Matched 

filtering  of  the  differenced  output  currents  then  yields  a  measurement  of 


M  1  a 


j0  [^(t)  -  12(t)]dt  ,  (85) 

where  normalization  by  the  electron  charge  has  been  used,  for  convenience,  to 
make  the  observation  values  dimensionless.  We  assume  that  the  mean  function 
and  covariance  function  of  the  stationary  complex-Gausslan  local-oscillator 
random  process 


lit)  i  (PL0(t)/hvoAd)1/2  exp(-Jew(t)) 


(86) 


are 


m  ■  C(l-Y)<NLQ>/AdT] 


1/2 


(87) 


and 
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K^(x)  -  (y<NL0>/AdT)k(T)  , 


(88) 


respectively,  where  <N^q>  Is  tha  average  number  of  local  oscillator  photons 

present  ovar  AdxT  ,  y^2  Is  tha  fractional  root-mean-square  (ms)  local 
oscillator  amplitude  fluctuation  level ,  and  k(r)  Is  a' real -valued 
normalized  covariance  (k(0)  ■  1).  It  than  follows  that 


<M>  ■  n(l-2e)<N1>g> 

+  4n[e(1-c)(l-Y)<NLg>]1'/2  <a?1>  , 


*sr 


(89) 


and 


<aM2>  ■  n(l-n)<N^g>  +  n2(1-2e)2<M^g>  +  16n2e(l-e)(l-Y)<N|_Q><4aji2> 


+  Cn(l-2e)<NL0>/T]2 


dT[Y2k2(r)  +  ZY(1-Y)k(T)]  (T-M) 


+  (2n/T)2(1.2c)YCc(l-e)<NL0>3(l-Y)]1/22<asi> 

T 


drk(r) (T- 1 r | ) 


+  (2n/T)2c(l-c)v<NL0>[T2+2<as+as> 


drk(T)(T-|r|)]  , 


(90) 


give  tha  mean  and  variance  of  the  M  measurement  from  which  a  signal-to-nolse 
ratio 

-  2 


SNRj^  a  <M>  /<amS 


(91) 


may  be  calculated.  In  Eq.  (89),  tha  first  tern  Is  tha  average  local  oscillator 
bias  current  contribution,  and  tha  second  term  Is  the  average  signal  field 
mixing  tern  contribution,  where  a^  ■  (a^+ag  )/2  for  a^  the  annihilation 
operator  of  tha  sola  excited  Fj  mode.  In  Eq.  (90),  tha  first  term  Is  due  to 
the  n<l  noises  nyacl  and  rvaC2  •  the  second  term  Is  due  to  the  local 
oscillator  quantum  noise  n^  ,  the  third  term  Is  the  signal  field  quantum 
noise,  tha  fourth  term  is  due  to  the  local -oscillator  power  fluctuations,  and 
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the  last  terms  are  due  to  the  random  modulation  of  the  mixing  current  by  the 
local  oscillator  fluctuations. 

The  previous  dual-datactor  homodyne  studies  of  Yuen  and  Chan  [25]  and 
Schumaker  '[30]  assume  c  *  1/2,  ( T -y)<N^q»>>1 »  and  a  slowly  fluctuating  local 
oscillator  (corresponding,  In  our  case,  to  k(t)  ■  1  for | t ) <T) .  In  this  limit  both 
prior  studies  find  (in  our  notation) 


(92) 


For  small  fractional  rms  local -oscillator  amplitude  fluctuations  (y«1), 

(93i  differs  from  (921  because  of  an  additional  noise  term  in  the  denominator 

.  +, 

that  Is  approximately  y(l  +  2<aj  a^> )/4 .  Physically,  this  term  arises  from  the 

random  local -oscillator  modulation  of  the  mixing  current,  an  effect  neglected  by 

the  earlier  studies.  In  order  for  this  term  to  be  Insignificant  compared  to  the 

»  2 

signal  quantum  noise  of  a  coherent  state  (<ca^  >  ■  1/4),  we  require  that 
•  +. 

v<as  as>  «1  ,  (94) 

1.*,,  the  fractional  rms  local  oscillator  amplitude  fluctuation  must  be  much 

smaller  than  the  square  root  of  the  reciprocal  of  the  average  number  of  signal 

field  photons!1  This  requirement  becomes  even  more  stringent  If  a  squeezed  state 

*  2 

Is  being  probed,  for  which  <aas1  >  <1/4  prevails. 


In  addition  to  exhibiting  the  potentially  significant  random  modulation 
of  the  mixing  term,  our  formulation,  (8**),  shows  another  effect  suppressed 
In  [25]  and  [30].  This  Is  the  random  amplitude  modulation  of  the  local  oscillator 
and  n  <  1  quantum  noises  by  the  classical  amplitude  noise  of  the  local 
oscillator.  Although  this  modulation  does  not  explicitly  enter  the  signal -to- 
nolse  ratio  ,  It  does  make  the  last  three  terms  In  (84)  non-Gausslan  random 
processes,  an  effectwhlch  will  modify  digital  communication  error  probability 
calculations  somewhat. 

Heterodyne  Detection  The  semi  cl  ass  leal  description  for  the  differenced  output 
currents  In  heterodyne  detection  Is 

1 1  (t)  -  12(t)  »  eri(l-2e)Pl0(t)/h\Jo 

+4en[c(l-c)PL0(t)/hvJoAd]1/2Re(j  dxFs(x.t)e‘J(2,rvIFt*,>L0(t))) 

Ad 

+e[n(l-c)PL0(t)/hvo]1/2  nJhotl(t) 

■«Cn«PLO(t)/hv0]%ot2(t)  ,  05) 

and  the  quantum  description  Is 

l^t)  -  12(t)  -  en(l-2c)PLQ(t)/hva 

Menlc(li)PLO(t)/hv0  Ad]1/2Re(f  dx^s(x,t)e"j(2ffVIFt‘*L0(t))) 

Ad 

+er,(l-2e)(PL0(t)/hvo)1/2nLQq(t) 

+e[n(l-n)(l-t)PL0(t)/hvo]VZnvacl(t) 
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-«Cn(1-n)cPL0(t)/hvo]1/2nvac2(t) 

+«ni;2cO-c)PL0(t)/hvo31/2nrq(t)  .  (96) 

Th««  results  differ  from  tha  corresponding  homodyne  results,  (83)  and  (84),  In 
only  two  respects.  First,  the  mixing  tarns  (second  tarms  on  tha  right  In  (95) 
and  (96))  baat  the  signal  flald  to  an  Intarmadlata  frequency  not  baseband, 
and  so  they  sense  both  quadratures  of  tha  signal  field.  Second,  the  quantum 
result  (96)  gains  a  noise  contribution  from  the  Image  band  quantum  noise  through 
n^(t)  .  The  local  oscillator  excess  noise  (and  Its  cancellation  when  c«  1/2) 
and  the  random  modulation  of  the  signal  and  noise  terms  by  the  local  oscillator 
fluctuations  thus  continue  to  be  present  In  the  heterodyne  case,  l.e.,  the 
Interpretations  given  for  tha  homodyne  situation  apply  here  as  well.  Once  again, 
the  relevant  previous  work  on  dual  detector  systems  [25],  [28],  [29]  does  not 
Include  all  the  effects  contained  In  our  treatment;  the  random  local  oscillator 
modulation  of  the  signal  and  noise  terms  Is  absent  In  the  above  analyses. 

As  an  Illustration  of  these  omissions,  let  us  compare  our  semi  classical  answer  (96) 
assuming  a  deterministic  monochromatic  plane-wave  pulse  signal  Fs(“,t)  ■ 

«S<W]/Z  for  xiA^.teT,  with  the  corresponding  t«l/2,  equal  quantum 

efficiency  result  of  Abbas  and  Chan  [29].  The  latter  claim,  In  our  notation,  that 
the  differenced  output  currents  consist  of  a  mean  current 

<(11(t)  -  12(t) )>-  2enT'1<NL0>1/2Re(asa’J2ffVIFt)  ,  (97) 

embedded  In  an  additive  zero-mean  white  Gaussian  noise  process  with  bilateral 
spectral  density 

S(f)  -  e2n<NLQ>/T  .  (98) 

We  have,  using  (86)  -  (87)  In  (95),  that  the  differenced  output  currents  consist 
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of  a  mean  current 

<(1-, (t)  -  12(t))>-2enT"1LO-Y)<NL0>31/2Re(ase'J2TTVIFt)  ,  (99) 

Plus  a  conditionally  non-statlonary  2ero-mean  white  Sausslan  shot  noise  process 
that,  gl yen  the  local  oscillator  power  waveform,  has  covariance  function 

K(t,s)  ■  e2n(PLQ(t)/hv0)$(t-s)  ,  (TOO) 

Plus  a  signal  dependent  zero-mean  stationary  Gaussian  noise  process 

Ht>  •  2en(Ad/T)1/2ReLc,s(l(t)  -  m^e'^IF*]  t  COT) 

with  covariance  function 

K1 (t )  ■  2(en/T)2v<NLQ>|as|2k(r)  cos(2it\>ift)  .  (102) 

When  Eq.(IOO)  Is  averaged  over  the  P^q  statistics  It  reduces  to  a  stationary 
white  noise  spectrum  (98)  ,  however  the  random  PLQ  fluctuations  make  the 
noise  non-Gausslan,  albeit  in  a  minor  way  If  y«l,  The  noise  current 
1"(t)  comas  from  the  random  modulation  of  the  mixing  term  and  may  present  a 
significant  degradation.  Consider  a  high  quality  (y«l),  slowly  fluctuating 
local  oscillator  (k(t)  5  1  for  |t|^T)  and  the  matched  filter  processor  generating 

M  "  8*1J0  1 2^ J2"1  /2cosC2tv  rfrt  -  arg(as)]  ,  (103) 

then  the  Abbas  and  Chan  model  gives  a  slgnal-to-nolse  ratio 

SNRm  ■  2n |as | 2  ,  (104) 


whereas  we  have  that 


2 


(105) 


snrm  - 


2n  lag 

1  +■  nvlaj 


As  In  the  quantum  homodyne  example  given  earlier,  at  high  average  detected 
signal  levels  there  Is  a  very  stringent  requirement  on  local  oscillator 
amplitude  fluctuations  If  SNR  degradation  Is  to  be  avoided, 

IV.  DISCUSSION 

At  this  point,  we  have  clearly  established  how  the  quantum  theory  for 
coherent  optical  detection  subsumes  the  familiar  semlclasslcal  statistics  In 
a  natural  way,  We  have  also  seen  that  the  quantum  approach  Is  essential  for 
studying  the  photodetection  statistics  of  non-classical  field  states.  There 
is  now  considerable  Interest  in  a  particular  class  of  non-classical  states, 
called  the  two-photon  coherent  states  [14]  or  the  squeezed  states  [15],  These 
states  are  In  essence  minimum  uncertainty  product  states  for  the  quadrature 
components  of  the  photon-units  field  operator  E(x,t).  In  particular,  for  a 
single  field  mode  with  annihilation  operator  a,  the  two-photon  coherent  state 
1 0 u , v  >  obeys  the  elgenket  relation 

(ua  +  vaT)|a;u,v>  ■  ajs;u,v>,  (106) 

where  a.u.v  are  complex  numbers  and  u,v  satisfy  |u|^  -  | v | ^  ■  1.  With 

a-|  ■  (a  +  a  +  )/2  and  a2  ■  (a  -  a+)/2j  denoting  the  quadrature  components  of  a, 

we  then  find  that  the  state  | a ; u , v >  gives 

<4«12>  “  (107a) 

and 

< Aa >  *  | u  +  v|^/4.  (107b) 
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When  u.v  a re  real  valued,  (107)  Implies  that  |s;u,v>  satisfies  the 
Heisenberg  relation 

<aa^>  *  &2  >  >  1/16  (108) 

with  equality,  as  does  the  familiar  coherent  state  |<*>  .  Unlike  the  coherent 
state,  which  gives  <aa^>  ■  <aa2^>  ■  1/4,  (107)  shows  that  there  Is  an 
asymmetric  noise  division  between  the  quadratures  (a  noise  squeezing)  In  the 

state  |a;u»v>  ,  with  the  low-noise  quadrature  being  less  noisy  than  a  coherent 
state.  This  noise  reduction  can  be  used,  In  principle,  to  effect  Important 
performance  improvements  In  optical  communications  [16]  -  [19]  and  precision 
measurements  [20]  -  [23], 

As  yet,  there  have  been  no  experimental  observations  of  squeezed  state 
light.  Theoretical  studies,  which  employ  varying  degrees  of  Idealization, 
indicate  that  such  states  may  be  generated  by  degenerate  four-wave  mixing  (DFWM) 
[39]  -  [42],  as  well  as  a  number  of  other  nonlinear  optical  processes  [14],  [15], 
[43]  -  [46].  We  are  presently  working  on  a  continuous -wave  DFWM  experiment  using 
homodyne  detection  to  generate  and  verify  the  quadrature  noise  squeezing.  In 
this  experiment,  a  single  frequency-stabilized  laser  will  be  used  to  provide 
all  the  Input  beams  to  the  four-wave  mixer,  as  well  as  the  local  oscillators 
for  dual -detector  homodyne  detection.  The  results  of  this  oaper  permit  the 
expected  photocurrent  statistics  for  this  experiment  to  be  derived,  Including 
the  effects  of  the  laser's  residual  amplitude  and  phase  fluctuations.  Specifi¬ 
cally,  an  Iterated  expectation  approach  is  used,  as  In  Section  III.  The  photo- 
current  statistics  are  first  obtained  assuming  the  laser  output  to  be  a  par¬ 
ticular  coherent  state.  This  entails  a  calculation  of  the  four-wave  mixer 
output  state,  along  the  lines  of  [40],  followed  by  a  calculation  of  the  sort 
performed  here  in  Section  IIIA.  To  average  over  the  input  laser  fluctuations , 
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we  assign  to  the  coherent  state  value  for  thl3  laser  a  classical  probability 
distribution.  We  can  then  proceed  as  In  Section  I I IB,  except  that  the  state 

of  the  signal  field  operator  In  the  homodyne  apparatus  Is  now  dependent  on  the 

/ 

coherent-state  value  of  the  local -oscillator  field  In  that  apparatus,  because 
both  fields  are  derived  from  the  same  laser. 
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Footnoter 

1.  The  convention  we  use  for  this  Fourier  transform  Is  necessitated  by 
the  accepted  quantum-optics  definition  for  what  constitutes  a  positive- 
frequency  field. 

2.  For  a  photomultiplier  tube,  the  Internal  current  gain  amplifies  the  current 

we  are  analysing  by  a  sufficient  amount  to  warrant  Its  treatment  as  a  classical 
entity.  In  the  coherent  optical  detection  cases  that  follow,  the  mixing  gain 
produced  by  the  strong  local  oscillator  has  a  similar  effect,  see  [22]. 

3.  A  classical  state  Is  either  a  Glauber  coherent  state  or  a  classically  random 
mixture  of  such  states.  In  either  ease,  the  density  operator  o  has  a 
proper  P-represantatlon  (17).  The  terminology  arises,  see  below,  because  a 
classical  state  p  gives. rise  to  the  same  statistics  In  quantum  photodetection 

‘  theory  as  found  for  a  classical  field  In  semi  classical  photodetectlon  theory, 

4.  Because  our  Idealized  detector  model  neglects  Internal  noise  sources  (dark 
current  shot  noise,  thermal  noise,  etc.)  N  from  Eq.  (El)  corresponds  to  the 
output  of  a  pulse-dlscrlmlnator/counter  applied  to  the  output  current 

j  dx  J(x,t).  In  other  words,  Eq.  (21)  models  the  output  of  an  Ideal  (unity 
Ad 

quantum  efficiency)  photomultlpllar-tuba/pulse-counter  setup. 

5.  Our  choice  for  the  beam  spllttor  transformation  agrees  with  that  employed  In 
[18],  and  Implies  that  the  field  leaving  the  other  port  of  this  optical  element 
Is  -(1-c)1/^Es(7,t)+e''/^SLQ{x,t) .  Other  beam  splitter  relations  (see,  e.g. 
[25]  ,  [30])  are  equivalent  to  ours  after  redefinition  of  the  Input  and 
output  planes. 

6.  A  critical  aspect  of  the  strong  local  oscillator  condition  acting  through  the 
measurement  operator  (48)  Is  that  the  mean  local  oscillator  field  and  Its 
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quantum  nolsa  both  contribute  to  Jh0ffl  through  the  direct  detection 
term,  but  only  the  mean  local  oscillator  field  (not  Us  quantum  noise) 
contributes  to  through  the  mixing  term. 

7.  Very  Interesting  noise  reductions  can  accrue  when  the  signal  and  image 
bands  are  quantum-mechanical ly  dependent  [22],. 

8.  Because  the  In-phase  and  quadrature  components  of  Fs  are  non-commuting 
observables,  the  Image  band  noise  enters  Into  heterodyning  In  order  to  enforce 
the  Heisenberg  uncertainty  principle  on  Ideal  simultaneous  observations  of 
these  Incompatible  quantities  (see  [36],  [37]), 

9.  Implicit  In  this  conditioning  statement  Is  the  fact  that  the  local  oscillator 
must,  with  very  high  probability,  remain  sufficiently  strong  to  ensure  the 
validity  of  the  Section  II  theory,  Also  note  that  the  signal  field  statistics 
may  depend  on  the  value  of  the  local  oscillator  field,  such  as  occurs  In  a 

.  laboratory  experiment  when  the  same  laser  Is  used  to  obtain  both  the  signal 
and  local  oscillator  beams  (see  Section  IV). 

10.  The  local  oscillator  fluctuations  must  not  be  such  as  to  Invalidate  the 

* 

Section  II  theory  for  any  state  | «j_q>  that  occurs  with  appreciable  probability. 
Also,  the  signal  state  (density  operator)  may  depend  on  the  value  of  the  local 
oscillator  field,  If,  for  example,  both  beams  originate  from  the  same  laser 
(see  Section  IV), 

11,  For  example,  to  keep  this  added  noise  below  10*  (In  standard  deviation)  of 

the  coherent-state  signal  quantum  noise  whan  <a^  a^>  *  10  ,  we  can  tolerate  no 

more  than  0.3%  local  oscillator  amplitude  fluctuation.  This  limitation  may 
be  significant,  in  precision  measurement  apollcatlons  for  which  signal -to-nolsa 
ratios  far  In  excess  of  40dB  are  sought. 
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Figure  Captions 

Fig.  1  Geometry  of  an  idealized  surface  photoemitter  of  active  region  Ad. 
Fig..  2  Configuration  for  optical  homodyne  detection. 

Fig.  3  Configuration  for  optical  heterodyne  detection. 

Fig.  4  Configuration  for  dual-detector  coherent  optical  detection. 
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